New class of g-modes and unexpected convection in neutron stars 
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We suggest a specific new class of low-frequency g-modes in superfluid neutron stars. We de- 
termine the Brunt- Vaisala frequency for these modes and demonstrate that they can be unstable 
with respect to convection. The criterion for the instability onset (analogue of the well known 
Schwarzschild criterion) is derived. It is very sensitive to equation of state and a model of nu- 
cleon superfluidity. In particular, convection may occur for both positive and negative temperature 
gradients. Our results have interesting implications for neutron star cooling and seismology. 
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Introduction. — This Letter is devoted to gravity os- 
cillation modes (g-modes) and the related phenomenon 
of convection in neutron stars (NSs). The restoring force 
for g-modes is buoyancy that originates from the depen- 
dence of the pressure on at least two quantities (e.g., 
density and temperature or density and chemical compo- 
sition), g-modes and convection are actively studied in 
laboratory experiments (e.g., [IrQl) and are widespread 
in nature. For instance, g-modes are observed in Earth 
atmosphere and ocean, in white dwarfs 4 \ , in slowly pul- 
sating B-stars [5|, and in other objects 6], while convec- 
tion is typical for most of the stars (including the Sun). 
In application to NSs, g-modes were studied, e.g., in Refs. 
|7Hl2j . In all these works the NS matter was assumed to 
be nonsuperfiuid (nor mal) . However, according to mi- 
croscopic calculations 0, 0], baryons in the internal 
layers of NSs become superfluid (SFL) at temperatures 
T < 10 8 -f- 10 10 K which has a drastic impact on stel- 
lar dynamics and evolution 0| . Recent real-time 
observations [16[ of the cooling NS in Cas A supernova 
remnant indicate that this NS has an SFL core |17l . 0] . 
A number of attempts U21 have been made to theo- 
retically predict g-modes in cold SFL NSs, but they have 
failed. This led to a general belief that g-modes do not 
exist in SFL interiors of NSs or, more precisely, their 
frequencies u are all degenerate at zero. In this Letter 
we show that proper account of finite temperature ef- 
fects extracts g-modes from the zero frequency domain. 
Moreover, these modes can be unstable with respect to 
convection. Possible applications of these results are out- 
lined. Below the Planck constant, the speed of light, and 
the Boltzmann constant equal unity, H — c — k-g, = 1 . 

Convection in NSs and SFL g-modes. — For simplic- 
ity, consider npe NS cores, composed of neutrons (n), 
protons (p), and electrons (e). To start with, assume 
that all particles are nonsuperfiuid. Any thermodynamic 
quantity in npe-matter (e.g., the heat function w — e+P, 
where e is the energy density and P is the pressure) 
can be presented as a function of 3 variables, say, P, 
x e = n e /nb, and xs = S/nj,. Here rii is the number den- 
sity for particles i = n, p, and e; n& is the baryon number- 
density; 5 is the entropy density. What is the local cri- 



terion for the absence of convection in npe-matter? It 
is easy to derive it in the same manner as it was done 
in Refs. 22, 0| (see also 11 1). Assume that a spheri- 



cally symmetric star is in hydrostatic equilibrium (but 
not necessarily in thermal or beta-equilibrium), that is 
VP = — w V0, where </>(r) is the gravitational potential 
and r is the radial coordinate. Here and below V = d/dr 
because all quantities of interest depend on r only. Con- 
sider two close points 1 and 2 with r = r\ and r 2 . Let A\ 
and A2 be the values of some thermodynamic quantity 
A at points 1 and 2, respectively, and AA = A% — A\. 
Let us displace adiabatically a small fluid element up- 
ward from point 1 to point 2. At point 2 the pres- 
sure of the fluid element adjusts itself to the surround- 
ing pressure Pi = P\ + AP, while x e and xs remain 
unchanged and equal to x e % and xsi (we assume that 
beta-processes are slow). The matter is stable against 
convection if the inertial mass density of the lifted ele- 
ment is larger than the equilibrium density at point 2. 
For relativistic matter the role of the inertial mass den- 
sity is played by w 0. Thus, the stability requires that 
w(Pi, x e2 , XS2) <w(P 2 , x e i, xsi). Expanding w in Tai- 
lor series near point 1, we obtain 

d Xc w(P, x e , Xs) Vi e + d Xs w(P, x e , xs) Vxs < 0, (1) 

where &a = d/dA. When w is a function of P and xs 
only, Eq. ([T]) immediately reproduces the textbook cri- 



terion for the absence of convection (see, e.g., |22j). In 
a strongly degenerate matter the second term in Eq. (|T]) 
is much smaller than the first one and can be neglected. 
Similarly, to calculate the first term in Eq. (JlJ it is suf- 
ficient to set T = and xs = 0. Then, Eq. ([1]) reduces 
to d Xe w(P, x e ) \7x e < 0. This Ledoux-type criterion is 
always satisfied in beta-equilibrated NSs, i.e., they are 
stable against convection. Oscillations of such a matter 
near equilibrium correspond to temperature-independent 
composition g-modes, first studied in Ref. [Ioj |. 

Assume now that neutrons (and possibly protons) are 
SFL. What will be the analogue of criterion ([1]) in that 
case? Nucleon SFL leads to the appearance of two in- 
dependent velocity fields: SFL neutron velocity V sn and 
velocity of normal liquid component V q , composed of 
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neutron Bogoliubov excitations, protons and electrons 
[25j . The presence of extra velocity field V sn results 
in additional (besides equation VP = — wS7<j> ) c ondi- 
tion of hydrostatic equilibrium in SFL matter 2(| 27 1: 
V ul«„ e'H = 0, where fi n is the relativistic neutron chem- 
ical potential. As a result, when we displace the fluid el- 
ement, 'attached' to the normal liquid component, from 
point 1 to point 2, both P and \i n adjust themselves 
to their equilibrium values P 2 and fi n2 at point 2. The 
pressure adjusts by contraction/expansion of the fluid el- 
ement, while fi n adjusts by the variation in the number of 
'SFL neutrons' in this element. Note that, since SFL neu- 
trons can freely escape from the fluid element attached 
to the normal particles, the total number of neutrons in 
the element is not conserved, and neither are the quan- 
tities x e — n e /rib and xs — S/n b - In this situation the 
conserved quantity is x e s — S/n e , because both the en- 
tropy a nd electrons flow with the same velocity V q (e.g., 
26j-|28[). Bearing this in mind, it is convenient to con- 
sider w as a function of P, and x e s- Then the con- 
dition for stability against convection can be written as 

w(P 2 , Mn2, XeS2) < w(P 2 , fl n2 , XeSl) Or 







>z cS w(P, Mn, Xes) Vx eS < 0. 



(2) 



A similar condition was derived in a different way in Ref. 
[29| , where internal gravity waves were analyzed in a mix- 
ture of SFL He-4 and a normal fluid (see also 3fl 3lj|). 
Note that the left-hand side of Eq. @ depends on T and 
vanishes at T — 0. Then the system is marginally sta- 
ble, since there is no restoring force acting on a displaced 
fluid element. Thus, it is not surprising that the authors 
of Refs. EiSl, who assumed T = 0, did not find g- 
modes in SFL NSs. In contrast, consistent treatment of 
the temperature effects should reveal g-modes. 

To check it we perform a local analysis of SFL hydro- 
dynamic equations (see, e.g., [2M [SS]), describing oscilla- 
tions of an NS in the weak- field limit (cj) <C 1) at T ^ 0. 
We analyze short-wave perturbations, proportional to 
exp(iwi) exp[j J dr'k(r')]Yi m , where the wave number 
k of a perturbation weakly depends our(fc> \d\nk/dr\, 
WKB approximation), and Y; m is a spherical harmonic. 
Solving oscillation equations in the Cowling approxima- 
tion (in which (j) is not perturbed [13]), we find the stan- 
dard short-wave dispersion relation for the SFL g- 
modes: u? = M 2 1(1 + l)/[l(l + 1 ) + k 2 r 2 }, where 



M* = — 



(i + y) 



y 



d XcS w(P, fl n , Xes) VXeS (3) 



is the corresponding Brunt- Vaisala frequency squared; 
g = V0; y = n b Y pp /[p n (Y nn Y pp -Y 2 p )]-l > OJ-^ being 
the relativistic entrainment matrix (sec, e.g., 33, 34] and 
comment [35j]). The stability condition for these g-modes, 
N 2 > 0, coincides with Eq. (0. 

For numerical evaluation of Eq. © it is convenient to 
introduce a new set of independent variables rib, n e , and 



T instead of P, \x n , and x e s- Then Eq. ([3]) is approxi- 
mately rewritten as 



JV 2 w a C V F [- VT/( 5 T) + F - 1] 



(4) 



where we neglect small terms of the second and higher 
orders in T/e-p (ep is the typical particle Fermi energy). 
In Eq. flU) a = g 2 T (1 + y)/(y (J, n n b ) > and C v = 
Td T S > 0. Finally, F = 1 + ix n G 2 /{Gx C v ), where 
Gi = <9„ e P d nb Hn ~ d n „P d ne fi n and G 2 = <9„ e P d nb S + 
dn^fJ-n (S — ribd nb S — n e dn e S). For an NS in thermal 
equilibrium, the red-shifted temperature T°° = Te^ is 
constant throughout the core, VP 00 = (VP+.gP) = 0. 
In that case M 2 in Eq. (j4]) is positive and reduces to 



M 2 w a C v F 2 



(5) 



If NS is not in thermal equilibrium, M 2 can be neg- 
ative for certain VT°°. These gradients follow from 
Eq. ((4]) (or Eq. [2]) and are defined by the inequality 
FVT 00 > gT^F 2 , which is the analo gue of the ordinary 
Schwarzschild criterion for convection [23] . This inequal- 
ity, as well as Eq. ©, is valid not only in the weak-field 
limit ^<1, but also in the full general relativity. When 
it is satisfied, convective instability occurs [36]. Thus, 
the critical gradient for the instability onset is given by 



VT^ t = 5 T°°P. 



(G) 



Clearly, F determines the sign of VT c ^ t . If F < one 
gets the instability while heating the matter from below; 
if F > the instability occurs when it is heated from 
above. Note that both signs of VT°° can be realized in 
cooling NSs [37 1. 

Results. - First, consider the limit T <C T cn , T cp 
(T c i is the critical temperature for particles i = n, p), 
in which the nucleon entropy is exponentially suppressed 
and the entropy density S is provided by electrons, 
S = S e = T(3ttV) 2 / 3 /3. In this limit both J\f and 
VP c ^ t are oc T. Fig. QJa) presents the Brunt- Vaisala fre- 
quency Af (see Eq. [S]) as a function of rib for npe-matter 
in thermodynamic equilibrium for 5 equations of state 
(EOSs) and T = 10 7 K. One sees that for any EOS N 
vanishes at a certain n b that corresponds to F = 0. Fig. 
QJb) shows the critical gradient VT c ^- t (see Eq. [6]) versus 
n b for the APR EOS [!§]. As expected, VT c ^ t changes 
sign when F — 0. The region of parameters, where con- 
vection occurs, is filled with gray. 

When T is not too low (0.1P cl < P < T ci , i = n 
and/or p), nucleonic contribution to S is non-negligible. 
The results in that case strongly differ from those ob- 
tained in the limit T <C T cn , T cp and are presented in 
Fig. H For illustration, we adopt the APR EOS and 
take T = 1.5 x 10 s K. Some realistic profiles of singlet 
proton T cp (rib) and triplet neutron T cn (nt>) critical tem- 
peratures, employed in our numerical calculations, are 
shown in Fig. [DJa). Figs. [2fb,c) demonstrate N(rib) and 
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FIG. 1: Panel (a): J\f versus rib f° r EOSs of Armani et al. 
HJ, APR [13, and PAL [1J. We adopt the model I of PAL 
family with three values of the compression modulus, 120, 
180 and 240 MeV. Panel (b): VT^ versus n b for APR EOS. 
The instability region is filled with gray. Vertical dot-dashed 
line corresponds to the crust-core interface. Both panels are 
plotted for T = 10 7 K. Here and in Fig. [Ug = 10 14 cm s -2 . 

VT c ^ t (nf,), respectively. Solid lines in Figs. EIb,c) are 
obtained for T cn (rib) and T cp (rib) from Fig. EI a ). Dot- 
dashed lines are plotted for T cn (rib) from Fig. E{a), but 
for T cp — > oo. Finally, dashed lines in Figs. EJb,c) cor- 
respond to the limit T cn , T cp —> oo of Fig. [TJ (note, how- 
ever, that Fig. Q] is plotted for different T). For sim- 
plicity, we assume that the entropy density Si of parti- 
cles i = n, p, and e depends only on rii and T, so that 
S = S n (n„,T) + S p (n p ,T) + S e (n e ,T). Then F, which 
enters Eqs. ([5|) and (j6|), can be rewritten as 

F = 1 + 77~7T~ [( 9 n e P ~ n n d ne (i n ) d nn S n 

-n e <9„ c /i„ d np S p - n e d nc ix n d n<s S e + S d ne Hn] ■ (7) 

Clearly, the behaviour of M and VT c ^ t results from in- 
terplay between various derivatives of Si . The key role is 
played by the nucleon derivatives d ni Si (i — n, p). Be- 
cause Si strongly depends on T C i , which is in turn a very 
strong function of nb (especially, on the slopes of T C i(rib), 
see Fig. Et), d ni Si can be very large and thus determine 
the dependence of M and VT c ^ t on rib- 
To illustrate this point let us compare dot-dashed and 
dashed curves in Fig. EIc). Dashed curve is plotted as- 
suming S n = S p = (T c „, T cp — > oo) while for dot- 
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n b (fm 3 ) 

FIG. 2: Panel (a): T cn and T cp versus rib- Panels (b) and (c): 
M and VT c ^ t versus n b for EOS APR and T = 1.5 X 10 s K [see 
the horizontal dotted line in panel (a)]. Solid lines in panels 
(b) and (c) are obtained for T cn (rib) and T cp (nj,) from panel 
(a); dot-dashed lines: T cn (rib) is from panel (a), T cp oo; 
dashed lines: T cn , T cp — > oo. The right vertical dot-dashed 
line indicates the boundary between the SFL and normal neu- 
tron matter (in the latter SFL g- modes are absent). The ver- 
tical dotted line shows a similar boundary for protons. Other 
notations are the same as in Fig. [T] 

dashed curve only S p — 0. The terms in F related to 
dn n S n are negative for rib < 0.2 fm -3 and positive for 
ni, <^ 0.2 fm~ 3 . This can be easily understood if one 
bears in mind that for densities of interest (i) G\ > and 
d Uc P — n n d ni! fi n > 0; (ii) S n decreases with increasing 
T c „; and (Hi) T cn reaches maximum at rib ~ 0.2 fm 
(Fig. Et)- As a result, VT c ^ t given by the dot-dashed 
curve vanishes at lower rib than for dashed curve (Fig. 

Et). 

Now let us consider the effects related to the pro- 
ton entropy density S p (solid line in Fig. Et)- When 
T <C T cp the solid and dot-dashed curves coincide, be- 
cause in that case S p is negligible. When T approaches 
T cp with growing rib, VT c ^ t rapidly increases because 
dT C p/drib is large and negative (i.e. d np S p > 0) and 
dn^^n < [see the corresponding term in Eq. (|7J)]. At 
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T = T cp (vertical dotted line in Fig. [2j, d np S p and 8tS p 
are discontinuous that results in discontinuities of M and 
VT^. At rib <L 0.26 fm -3 all protons are normal and 
S p = T m*ppp/3, where m* andpFp are the proton effec- 
tive mass and Fermi momentum, respectively. At such 
rib the proton contribution to F is negative, i.e. the solid 
curve in Fig. [2jc) goes lower than the dot-dashed curve. 
The most important conclusion drawn from the analysis 
of Fig. [2] is that SFL g-modes and convection in the in- 
ternal layers of NSs are extremely sensitive to the EOS 
and the model of nucleon SFL. An account for the sin- 
glet neutron SFL at lower densities rib < 0.08 fm -3 may 
additionally affect M and VT c ^? t near the crust-core in- 
terface. 

Discussion and conclusion. - Our results indicate 
that NSs can have convective internal layers. This could 
affect the thermal evolution of young NSs (such as in Cas 
A), for which VT°° is not completely smoothed out by 
the thermal conductivity, as well as the thermal relax- 
ation of quasi-persistent X-ray transients 4l|, |42J . Note 



that in this Letter we only consider SFL g-modes and 
convection in the NS cores (but not in the crust). If the 
NS crust is elastic, as it is usually assumed, then it is 
most likely that core SFL g-modes do not penetrate the 
crust, while the crustal SFL g-modes are 'mixed' with the 
shear modes @ (for which the restoring force is elastic- 
ity) , and pushed to frequencies v ^> M (but see comment 
[431]). In that case convection is absent. However, if the 
inner crust (especially, mantle [Hj]) is plastic [H| then the 
existence of crustal SFL g-modes and convection cannot 
be excluded, that can have even more interesting impli- 
cations for NS cooling. 

We have shown that g-modes can propagate in SFL NS 
matter. But how can they be excited? Among the po- 
tential scenarios is the excitation of stable SFL g-modes 
by unstable ones (i.e., by convective motions). Another 
possibility was considered in Refs. [Ill El in applica- 
tion to composition g-modes of normal NSs. It consists 
in resonant excitations of SFL g-modes by tidal inter- 
action in coalescing binary systems, when the frequency 
of the tidal driving force equals one of the g-mode fre- 
quencies. Finally, SFL g-modes in rotating NSs could 
be excited due to gravitational driven (CFS) instability, 
though this scenario does not seem very realistic [l2| be- 
cause J\f is low, which results in a large gravitational 
radiation timescale. 

To conclude, we have predicted a new class of g-modes 
in SFL NSs. We have calculated their Brunt- Vaisala fre- 
quency A/", which strongly depends on T and vanishes 
at T = 0. The SFL g-modes appear to be unstable 
for certain temperature gradients (that correspond to 
Af 2 < 0). We have derived the criterion for convective 
instability (analogue of the Schwarzschild criterion) in 
SFL NS cores. We have shown that convection in the NS 
core may occur for both positive and negative temper- 
ature gradients and is extremely sensitive to the model 



EOS and nucleon SFL. We have only outlined the prop- 
erties of SFL g-modes. In particular, we have not calcu- 
lated their frequency spectrum and damping times. We 
plan to fill these gaps in the future publication. Though 
we have only considered npe-matter of NSs, our anal- 
ysis can be easily extended to SFL hyperon and quark 
stars, for which we also predict the existence of global, 
low-frequency SFL g-modes. 
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